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Abstract 

For the iron-based high T c superconductor LaFeAsOi_ x F x , we construct a minimal 
model, where all of the five Fe d bands turn out to be involved. We then investi- 
gate the origin of superconductivity with a five-band random-phase approximation 
by solving the Eliashberg equation. We conclude that the spin fluctuation modes 
arising from the nesting between the disconnected Fermi pockets realise, basically, 
an extended s-wave pairing, where the gap changes sign across the nesting vector. 
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1 Introduction 



While the physics of high-T c cuprate has matured after the two decades since 
the discovery, the superconductivity the iron-based pnictide LaFeAsO doped 
with fluorine discovered by Hosono's group pQ is more than welcome as a fresh 
challenge for yet another class of high-T c systems. Indeed, the iron-based ma- 
terial, along with various other ones in the same family of compounds with 
higher transition temperatures (T c )[2j, are remarkable as the first non-copper 
compound that has T c 's exceeding 50 K. This immediately stimulates renewed 
interests in the electronic mechanism of high T c superconductivity. In order to 
investigate the pairing mechanism, here we first construct an electronic model 
for LaFeAsOi-^Fu; using maximally localised Wannier orbitals obtained from 
first principles calculation. The minimal model turns out to involve all the five 
Fe d orbitals. [3] Hence the iron-based material is contrasted with the cuprate, 
which is a one-band, doped Mott insulator. We then apply the random-phase 
approximation (RPA) to solve the Eliashberg equation. We conclude that a 
nesting between multiple Fermi surface (pockets) results in a development of 
a peculiar spin fluctuation mode, which in turn realises an unconventional 
pairing, which is basically an extended s-wave where the gap function changes 
sign across the nesting vector. [3P] The result is intriguing as a realisation of 
the general idea that the way in which electron correlation effects appear is 
very sensitive to the underlying band structure and the shape of the Fermi 
surface. [5] 
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2 Band structure 



LaFeAsO has a layered structure, where Fe atoms form a square lattice in each 
layer, which is sandwiched by As atoms. [6] Due to the tetrahedral coordina- 
tion of As atoms, there are two Fe atoms per unit cell. The experimentally 
determined lattice constants are a = 4.03A and c = 8.74A, with two internal 
coordinates = 0.142 and za s = 0.651. [7] We have first obtained the band 
structure (FigJTJ) for these coordinates with the density-functional approxima- 
tion with plane- wave basis [8], which is then used to construct the maximally 
localised Wannier functions (MLWFs) [9 J . These MLWFs, centered at the two 
Fe sites in the unit cell, have five orbital symmetries [d^-W--, d X z, dyz, 
d x i-Y 2 -, dxYi where X, Y, Z refer to those for this unit cell with two Fe sites 
as shown in the bottom panel of FigJT]). The two Wannier orbitals in each unit 
cell are equivalent in that each Fe atom has the same local arrangement of 
other atoms. We can then take a unit cell that contains only one orbital per 
symmetry by unfolding the Brillouin zone. [TO] and we end up with an effective 
five-band model on a square lattice, where x and y axes are rotated by 45 de- 
grees from X-Y. We refer all the wave vectors in the unfolded Brillouin zone 
hereafter. We define the band filling n as the number of electrons/number of 
sites (e.g., n = 10 for full filling). The doping level x in LaFeAsOi^F,!. is 
related to the band filling as n = 6 + x. 

The five bands are heavily entangled as shown in Figj2] (left panel) reflecting 
strong hybridisation of the five 3d orbitals, which is physically due to the 
tetrahedral coordination of As atoms around Fe. Hence we conclude that the 
minimal electronic model requires all the five bands. In Figj2](right), the Fermi 
surface for n = 6.1 (corresponding to x — 0.1) obtained by ignoring the inter- 
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layer hoppings is shown in the two-dimensional unfolded Brillouin zone. 

The Fermi surface consists of four pieces (pockets in 2D): two concentric hole 
pockets (denoted here as a%, 02) centered around (k x , k y ) = (0, 0), two electron 
pockets around (vr, 0) or (0,7r) (P2), respectively. a« (/%) corresponds to 
the Fermi surface around the TZ line (MA in the original Brillouin zone) in the 
first-principles band calculation. [TT] Besides these pieces of the Fermi surface, 
there is a portion of the band near (jr, tt) that is flat and touches the Ep at 
n = 6.1, so that the portion acts as a "quasi Fermi surface (7)" around (n, n), 
which has in fact an important contribution to the spin susceptibility. As for 
the orbital character, a and portions of j3 near Brillouin zone edge have mainly 
dxz and dyz character, while the portions of /3 away from the Brillouin zone 
edge and 7 have mainly dx 2 ~Y 2 orbital character (Figj3l bottom panels). 

An interesting feature in the band structure is the presence of Dirac cones, i.e., 
places where the upper and the lower bands make a conical contact. [T2"irT3] 
The ones closest to the Fermi level lies at positions where the dx 2 -Y 2 and the 
dxz /dyz bands cross, just below the (3 Fermi surface. 

3 Many-body Hamiltonian and 5-band RPA 

We consider a two-dimensional model where the inter-layer hoppings are ne- 
glected. For the many body part of the Hamiltonian, we consider the standard 
interaction terms that comprise the intra-orbital Coulomb U, the inter-orbital 
Coulomb U', the Hund's coupling J and the pair-hopping J'. The many body 
Hamiltonian then reads 
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Fig. 1. Top panel: The band structure of in the original (folded) Brillouin zone, where 
the symbols represent the LDA result while curves the ten-band model derived with 
the maximally localised Wannier orbitals. Bottom panel: The original (dashed) and 
the reduced (solid) unit cells in real space with • (Fe), V (As below the Fe plane) 
and A (above). 
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Fig. 2. Left: The band structure of the five-band model in the unfolded Brillouin 
zone. Orbital characters of the bands are also indicated. The short solid arrow 
denotes the position of the Dirac cone. Right: Fermi surface for n = 6.1 in the 
unfolded Brillouin zone, where the original Bz is indicated by dashed lines. The 
blue (red) portions have strong X 2 — Y 2 (Y Z, ZX) character. Dashed blue curves 
represent the portions where a X 2 — Y 2 charactered band lies very close to the Fermi 
energy for n = 6.1, which actually become a Fermi surface for smaller doping. 
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where i,j denote the sites and n, v the orbitals, and tf^ is the transfer energy 
obtained in the previous section. The orbitals d 3 z 2 -R 2 , dxz, dyz, d x ^-Y 2 i 
dxY are labeled as v = 1,2,3,4, and 5, respectively. As for the electron- 
electron interactions, there have been some theoretical studies that estimate 
the parameter values. Some give U 3> J. [T4"irPo] while others U ~ J.JTB] We 



6 



assume here that U ^> J, and take the values U = 1.2, U' = 0.9, J = J' = 
0.15 throughout the study. These values are smaller than the values obtained 
in ref. fT^lfTB] because the self energy correction is not taken into account in 
the present calculation, so that small values of interaction parameters are 
necesssary to avoid magnetic ordering at high temperatures. 

Having constructed the model, we move on to the RPA calculation, where 
the modification of the band structure due to the self-energy correction is not 
taken into account. Multiorbital RPA is described in e.g. ref. |17|[T%] . In the 
present case, Green's function Gi m (k) (k = (k,iu n )) is a 5 x 5 matrix. The 
irreducible susceptibility matrix 

xl lM (q) = T,G h i 3 (k + q)G kl2 (k) (2) 



(li = 1, ...,5) has 25 x 25 components, and the spin and the charge (orbital) 
susceptibility matrices are obtained from matrix equations, 



Xs(q) 
Xc(q) 



X°(q) 



1 - S X °(q) 

X°(q) 
l + C X °(q) 



(3) 
(4) 



where 



U, U l 1= l 2 =l 3 = l 4 



U', -U' + J l 1 = l 3 ^l 2 = l 4 



J, 2U'-J, l 1 = l 2 ^l 3 = l 4 



J', J' h — h ^ h — h 



(5) 
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We denote the largest eigenvalue of the spin susceptibility matrix for iuo n = 
as Xs{k). The Green's function and the effective singlet pairing interaction, 

V s (q) = ls X s (q)S - \c X c {q)C + \{S + C), (6) 

are plugged into the linearised Eliashberg equation, 
T 

Whh(k) = -— J2 22 v hi2hh(q)Gi 2 i 5 { k - q)<Pi 5 k(k - q)G hk {q - k). (7) 

1 hhhk 

The 5x5 matrix gap function <pi m in the orbital representation along with the 
associated eigenvalue A are obtained by solving this equation. The gap func- 
tion can be transformed into the band representation with a unitary transfor- 
mation. The temperature is fixed at T = 0.02eV throughout the study, and 
32 x 32 Appoint meshes and 1024 Matsubara frequencies are taken. We find 
that the spin fluctuations dominate in magnitude over orbital fluctuations as 
far as U > U', so we can characterise the system with the spin susceptibility. 

4 Result: spin structure 

Let us first look at the key ingredient: the spin susceptibility in the top panel 
of FigJHl The susceptibility Xs has peaks around (k x ,k y ) = (vr,0), (0, 7r). In 
addition we note that there is a ridge extending from (it,tt/2) to (tt/2, 7r) 
around (71", 7r). To explore the origin of these spin structures, we show Xs3333 
and Xs4444 m the middle of FigJHl which represent the spin correlation within 
dyz and dx^-Y 2 orbitals, respectively. x«3333 peaks solely around (ir, 0) and 
(0, 7r), which reflects the nesting between the XZ, FZ-charactered portions of 
a and j3 Fermi pockets as shown in a bottom panel of Figj3l On the other hand, 
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Fig. 3. Top panel: Five-band RPA result for the spin susceptibility \s for U = 1.2, 
U< = 0.9, J = J' = 0.15, n = 6.1 and T = 0.02 (in eV). Middle panels: Diagonal 
components, Xs3333 e ft) and Xs4444 (right)) of the spin susceptibility matrix in the 
orbital representation (3 : YZ, 4 : X 2 — Y 2 ). Bottom panel: Nesting is shown for the 
Fermi surface for orbitals XZ, YZ (left) and for X 2 —Y 2 (right). Here the thickness 
of the Fermi surface represents the strength of the respective orbital character. The 
gray areas around the corners in the left panel indicates the 7 "quasi" Fermi surface. 
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Xs4444 has peaks around (ir, 0), (0, ir) and around (tt, 7r/2)/(7r/2, 7r) as well. The 
former is due to the nesting between the 7 quasi Fermi surface and the dx 2 -Y 2 
portion of the j3 Fermi surface as was first pointed out in ref. [1], while the lat- 
ter originates from the nesting between the dx^-Y 2 portion of the f3\ and $2 
Fermi surfaces. [3] The (ir, 0), (0, ir) feature is consistent with the stripe (i.e., 
collinear) antiferromagnetic order for the undoped case, which was suggested 
by transport and optical reflectance. |19j and further confirmed by neutron 
scattering experiments. [7] The stabilization of such an antiferromagnetic or- 
dering has also been pointed out in first principles calculations. |12|[T$1I2T)] 



5 Result: superconductivity 

The presence of multiple set of nesting vectors revealed above provides an 
interesting case of the gap function in a spin-fluctuation mediated supercon- 
ductivity, since multiple nestings can not only cooperate but also compete 
with each other. Namely, the a-(3 and 7-/? nestings tend to favour the sign 
reversing s-wave pairing, in which the gap changes sign between a and j3 with 
a fixed sign (i.e., full gap) on each of the Fermi surface. [1] On the other hand, 
(3\ —fo nesting tends to change the sign of the gap between these two pockets, 
which can result in either <i-wave pairing or an s-wave pairing with nodes on 
the j3 Fermi surface [2TH22] . For the band structure of LaFeAsO (obtained by 
using the experimentally determined lattice structure), the (71, 0) spin fluctu- 
ation dominates, and the sign-reversing s-wave with no nodes on the Fermi 
pockets dominates for the present set of parameter values. [22] In FigJH we 
show the gap function in the band representation for the third and the forth 
bands, which produce the 02 and (3 Fermi surfaces, respectively. A number of 
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theoretical studies that adopt effective two band models, [2311241125] the present 
five band model. [26|27II28] or a 16 band dp model[29] have also found that this 
sign reversing s-wave is a good candidate of the gap function in this material. 
The sign change in the s-wave gap is analogous to those in models studied 
by Bulut et ai, [30] and also by the present author for the disconnected Fermi 
surfaces [3"Tf3"2"] . It is also reminiscent of the unconventional s-wave pairing 
mechanism for Na x Co02 • yH 2 0[33] proposed by Kuroki et al. [31] 

We have to realise, however, that the gap can vary significantly along the 
f3 Fermi surface, and also between different pieces of the Fermi surface due 
to the multiorbital character of the system. In the present case, the gap for 
the <ix 2 -Y 2 -charactered portion of the (3 Fermi surface is about twice as large 
as that for the XZ, YZ-charactered portions, namely near the Brillouin zone 
edge of the f3 Fermi surface and on the a%, «2 Fermi surface. However, we find, 
in further calculations, that this gap variance is not universal, and depends 
strongly on the electron density and also on the details of the band structure. 
This is because the (3 — 7 nesting and thus the (it, 0) vs. (it, tt/2) spin fluctu- 
ation competition is sensitively affected by the position of the dx 2 -Y 2 portion 
of the band near (jr, tt) with respect to the Fermi level. The details of this 
band filling and band structure dependences will be published elsewhere. 

6 Conclusion 

To summarise, we have constructed a five-band electronic model for LaFeAsOi_ x F 
which we consider to be the minimum microscopic model for the iron-based 
superconductor. Applying a five-band RPA to this model, we have found that 
spin-fluctuation modes around (tt, 0), (0, tt) develop due to the nesting between 
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Fig. 4. Gap function in the band representation for band 2 (top) band 3 (middle) 
and band 4 for the same set of parameter values as in FigJSJ 
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disconnected Fermi surfaces. Based on the linearised Eliashberg equation, we 
have concluded that multiple spin fluctuation modes realises unconventional, 
extended s-wave pairing, where the gap changes sign across the nesting vec- 
tors. 

So the general picture obtained here is that the iron compound is a multi- 
band system having electron and hole pockets, as sharply opposed to the 
cuprate which is a one-band and nearly half-filled system with a simply con- 
nected Fermi surface. This poses a challenging future problem of elaborating 
respective pros and cons for the iron compound and the cuprate for supercon- 
ductivity. 
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